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Abstract 

Reliable calculations of financial risk require that the fat-tailed nature of prices 
changes is included in risk measures. To this end, a non-Gaussian approach to 
financial risk management is presented, modeling the power-law tails of the returns 
distribution in terms of a Student-i distribution. Non-Gaussian closed-form solutions 
for Value-at-Risk and Expected Shortfall are obtained and standard formulae known 
in the literature under the normality assumption are recovered as a special case. 
The implications of the approach for risk management are demonstrated through 
an empirical analysis of financial time series from the Italian stock market and in 
comparison with the results of the most widely used procedures of quantitative 
finance. Particular attention is paid to quantify the size of the errors affecting the 
market risk measures obtained according to different methodologies, by employing 
a bootstrap technique. 
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1 Introduction 



A topic of increasing importance in modern economy and society is the devel- 
opment of reliable methods of measuring and controlling financial risks. Ac- 
cording to the new capital adequacy framework, commonly known as Basel II 
accord [1], any financial institution has to meet stringent capital requirements 
in order to cover the various sources of risk that they incur as a result of 
their normal operation. Basically, three different categories of risk are of in- 
terest: credit risk, operational risk and market risk. In particular, market risk 
concerns the hazard of losing money due to the fluctuations of the prices of 
those instruments entering a financial portfolio and is, therefore, particularly 
important for financial risk management. 

In the financial industry today, the most widely used measure to manage 
market risk is Value-at-Risk (VaR) [2,3]. In short, VaR refers to the maximum 
potential loss over a given period at a certain confidence level and can be 
used to measure the risk of individual assets and portfolios of assets as well. 
Because of its conceptual simplicity, VaR has become a standard component 
in the methodology of academics and financial practitioners. Moreover, VaR 
generally provides a reasonably accurate estimate of risk at a reasonable com- 
putational time. Still, as discussed in the literature [3,4], VaR suffers from 
some inconsistencies: first, it can violate the sub-additivity rule for portfolio 
risk, which is a required property for any consistent measure of risk, and, sec- 
ondly, it doesn't quantify the typical loss incurred when the risk threshold is 
exceeded. To overcome the drawbacks of VaR, the Expected Shortfall (or Con- 
ditional VaR) is introduced, and sometimes used in financial risk management, 
as a more coherent measure of risk. 

Three main approaches are known in the literature and used in practice for 
calculating VaR and Expected Shortfall. The first method consists in assum- 
ing some probability distribution function for price changes and calculating 
the risk measures as closed-form solutions. This approach is called parametric 
or analytical and is easy to implement since analytical expressions can often 
be obtained. The parametric approach usually relies on the (log)normality as- 
sumption for the returns distribution, although some analytical results using 
non-Gaussian functional forms are available in the literature [5,6]. Actually, it 
is well known that empirical price returns, especially in the limit of high fre- 
quency, do not follow the Gaussian paradigm and are characterized by heavier 
tails and a higher peak than a normal distribution. In order to capture the lep- 
tokurtic (fat-tailed) nature of price returns, the historical simulation approach 
is often used as an alternative to the parametric method. It employs recent 
historical data and risk measures are derived from the percentiles of the dis- 
tribution of real data. This method is potentially the most accurate because 
it accounts for the real statistics of price changes but it is computationally 
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quite demanding (especially when applied to large portfolios) and absolutely 
depending on the past history of empirical data. A third approach consists 
in Monte Carlo simulations of the stochastic dynamics of a given model for 
stock price returns and in calculating risk measures according to Monte Carlo 
statistics. This method, however, requires very intensive simulations to achieve 
risk measures predictions with acceptable numerical errors. 

As a result of the present situation, reliable and possibly fast methods to calcu- 
late financial risk are strongly demanded. Inspired by this motivation, the aim 
of this paper is to present a non-Gaussian approach to market risk manage- 
ment and to describe its potentials, as well as limitations, in comparison with 
standard procedures used in financial analysis. To capture the excess of kurto- 
sis of empirical data with respect to the normal distribution, the statistics of 
price changes is modeled in terms of a Student-t distribution, which is known 
to approximate with good accuracy the distribution derived from market data 
at a given time horizon [3,7] and is widely used in the financial literature. In 
the econophysics literature, the Student-t distribution is also known as Tsal- 
lis distribution, emerging within the framework of statistical physics [8]. It 
has been shown in various studies [8,9] that the distribution of returns can 
be modeled quite well by a Tsallis distribution, which, for this reason, has 
been already used in a number of financial applications, ranging from option 
pricing [10] to risk analysis [11]. However, with respect to the investigation of 
Ref. [11], we include in our analysis the study of the Expected Shortfall and 
we present, in the spirit of a parametric approach, analytical expressions for 
the risk measures in order to provide accessible results for a simple practical 
implementation. At a variance of the recent calculation in Ref. [6] , where ana- 
lytical results for risk measures using Student-t distributions are presented, we 
critically investigate the implications of our non-Gaussian analytical solutions 
on the basis of an empirical analysis of financial data and we perform detailed 
comparisons with the results of widely used procedures. 

The paper is organized as follows. In Section 2 non-Gaussian closed-form ex- 
pressions for VaR and Expected Shortfall are derived as generalizations of the 
analytical formulae known in the literature under the normality assumption. 
It is also shown how the standard Gaussian formulae of the parametric ap- 
proach are recovered, in the appropriate limit, as a special case. In Section 3 
an empirical analysis of daily returns series from the Italian stock market is 
performed, in order to constrain the Student-t parameters entering the for- 
mulae of Section 2 and to describe the ingredients needed for the forecoming 
risk analysis. The latter is carried out in Section 4. The implications of the 
parametric non- Gaussian approach for VaR and Expected Shortfall are shown 
in Section 4 and compared with the results of the parametric normal method, 
of its improved version known as RiskMetrics methodology and of the histor- 
ical simulation. Particular attention is paid to quantify the size of the errors 
affecting the various risk measures, by employing a bootstrap technique. Con- 
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elusions and possible perspectives are drawn in Section 5. 



2 Non-Gaussian closed-form expressions for risk measures 

Value-at-Risk is referred to the probability of extreme losses in a portfolio 
value due to adverse market movements. In particular, for a given significance 
level V* (typically 1% or 5%), VaR, usually denoted as A*, is defined as the 
maximum potential loss over a fixed time horizon At. In terms of price changes 
AS, or, equivalenty, of returns R = AS/ S, VaR can be computed as follows 

V* = f A dAS P At (AS) = S f A /S dR P At (R), (1) 

J—oo J —CO 

where P At (AS) and P A t(R) are the probability density functions (pdfs) for 
price changes and for returns over a time horizon At, respectively. For financial 
analysts, VaR has become the standard measure used to quantify market risk 
because it has the great advantage to aggregate several risk component into 
a single number. In spite of its conceptual simplicity, VaR shows two main 
drawbacks: it is not necessary subadditive and it does not quantify the size of 
the potential loss when the threshold A* is exceeded. 

A quantity that does not suffer of these disadvantages is the so called Expected 
Shortfall (ES) or Conditional VaR (CVaR), E\ It is defined as 

E* = ± dAS (-AS) P At (AS) = A j~^' S dR (-R) P At (R), (2) 
with V* and A* as in Eq. (1). 

The standard approach in the financial literature [2,12] is to assume the returns 
as normally distributed, with mean m and variance a 2 , i.e. R ~ Af(m, a 2 ). In 
that case, VaR and ES analytical expressions reduce to the following closed- 
form formulae 

A* = -mS + aS V2 erfrT 1 ^*) (3) 

and 

E* = -mSo + exp{-[erfc- 1 (2^)] 2 }, (4) 

where So is the spot price and erfc" 1 is the inverse of the complementary error 
function [13]. 

However, it is well known from several studies, especially in the econophysics 
literature [3,7,14], that the normality hypothesis is often inadeguate for daily 
returns and, more generally, for high-frequency stock price variations. A better 
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Fig. 1. Convergence of the VaR (left) and ES (right) Student-i formulae toward 
Gaussian results when approaching the limit v — ► +00. 
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agreement with data is obtained using leptokurtic distributions, such as trun- 
cated Levy distributions or Student-t ones. Despite this interesting feature, 
the former family has the main disadvantage that it is defined only throught 
its characteristic function and we have no analytic expression for the pdfs [15]. 
Moreover, in order to compute the cumulative density function (cdf), which 
is a necessary ingredient of our analysis, we have to resort to numerical ap- 
proximations. For the reasons above, to model the returns, we make use of a 
Student-t distribution defined as 



B(u/2,l/2) [ a 2 + (jR _ m)2 ]^ 



i 



(5) 



where v G (1, +oo) is the tail index and B(v/2, 1/2) is the beta function. It is 
easy to verify that, for v > 2, the variance is given by a 2 = a 2 /(u — 2), while, 
for v > 4, the excess kurtosis reduces to k = 6/(u — A). Under this assumption, 
we obtain closed-form generalized expression for VaR and ES given by 



A* = -mS + aSo^/u - 2 



1 - A* 



A* 



and 



E* = -mS + 



a S [A *^ 



V*B(i>/2, 1/2) v — 1 

where A* ' '' ' " ' '~ J 

function, according to the definition of Ref. [13]. 



(6) 



(7) 



(2V*) and ^ is the inverse of the incomplete beta 



As shown in Fig. 1, we have checked numerically the convergence of formu- 
lae (6) and (7) to the Gaussian results (3) and (4), in the appropriate limit 
v — > +oo. We chose v = 2.75,3.5,4.5,100 and m — 0, crS = 1, but we 
checked that the value of these last parameters does not affect the conver- 
gence, as expected. As can be seen, the points corresponding to v — 100 are 
almost coincident with the Gaussian predictions, demonstrating that our re- 
sults correctly recover the Gaussian formulae as a special case. It is also worth 
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noting that each line, corresponding to a fixed u, crosses over the Gaussian 
one for a certain V*. Analogously, for a fixed V*, there exists a Across value 
whose line crosses the Gaussian result at that significance level. In the light of 
this observation, we report in Table 1 the values of f cr0 ss corresponding to a 
given V* for both VaR and ES. As can be observed, the growth of u cross with 
V* is very rapid for VaR, while for ES and for usually adopted significance 
values, f cr oss keeps in the interval [2.09,2.51]. From this point of view, VaR 
and ES are quite different measures of risk, since the crossover values for the 
latter are much more stable than those associated to the first one. This result 
can be interpreted as a consequence of ES as a more coherent risk measure 
than VaR. 
Table 1 

Values of v crossover for VaR and ES corresponding to different significance levels 
V*. 





1% 


2% 


3% 


4% 


5% 


Across (VaR) 
Across (ES) 


2.44 
2.09 


3.21 
2.18 


5.28 
2.28 


32.38 
2.38 


> 100 
2.51 



3 Empirical analysis of financial data 



The data sets used in our analysis consist of four financial time series, com- 
posed of V = 1000 daily returns, from the Italian stock market. Two series 
are collections of data from the Italian assets Autostrade SpA and Telecom 
Italia (from May 15 t/l 2001 to May 5 th 2005), while the other two correspond 
to the financial indexes Mib30 and Mibtel (from March 27 th 2002 to March 
13 th 2006). The data have been freely downloaded from Yahoo Finance [16]. 

Figure 2 shows a comparison between the historical complementary cumula- 
tive density function P> of the negative daily returns and two theoretical fits 
obtained using Gaussian and Student-t distributions. The parameters values 
of the fitted curves, as obtained according to the likelihood procedure de- 
scribed below, are displayed in Table 2. In principle, we could perform the fit 
according to different methods, but we have to balance between accuracy and 
computational time. Therefore, we estimate mean and variance as empirical 
moments, i.e. 

■y N-l 

m = NT, Rt-i (8) 

i=0 

and 

-I N-l 



a 



2 - * Y.( R t-i-™)\ (9) 
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where R = (Rt, . . . , Rt-N+i) is the A^-dimensional vector of returns. Using 
the above m and a values, we derive a standardized vector (with zero mean 
and unit variance) r = (r t , . . . ,r t -N+i), where r t -i = {R t -i — m)ja for i = 
0, . . . , iV — 1. In order to find the best value for the tail parameter v, we look 
for the argument that minimizes the negative log-likelihood, according to the 
formula 

N-l 



v = argmm v>2 



i=0 



(10) 



where the constraint v > 2 prevents the variance to be divergent and ls 
as in Eq. (5), with m = and a = \Jv — 2. This apparently simple optimization 
problem can not be solved analytically. In fact, the normalization factor in the 
Eq. (5) does depend on the tail index v in a non trivial way. Actually, the beta 
function B(vj1, 1/2) only admits an integral representation and therefore we 
implemented a numerical algorithm to search for the minimum. 

As shown in Section 2, the excess kurtosis k depends only on v and this pro- 
vides an alternative and more efficient way to estimate the tail parameter [11]. 
However, this approach forces v to be bigger than 4, while from Table 2 it can 
be seen that all the exponents obtained in the likelihood-based approach are 
smaller than 3.5. For this reason, the implementation of the excess kurtosis 
method is inadequate for the time series under study here. In order to test the 
robustness of our results, we also performed a more general three-dimensional 
minimization procedure over the free parameters (m,a,v). The multidimen- 
sional optimization problem was solved by using the MINUIT program from 
CERN library [17]. The obtained numerical results are in full agreement with 
the previous ones, but the process is more computationally burden and more 
cumbersone, since it requires a lot care in avoiding troubles related to the 
appearing of local minima in the minimization strategy. 

In Fig. 2 we show the cumulative distribution P > obtained using the empirical 
parameters of Table 2. As expected, we measure daily volatilities of the order 
of 1% and quite negligible means (~ 0.01%). The tail parameters fall in the 
range (2.9,3.5), thus confirming the strong leptokurtic nature of the returns 
distributions, both for single assets and market indexes. The quality of our fit 
clearly emerges from Fig. 2, where one can see a very good agreement between 
Student-t and historical complementary cdfs, while the Gaussian distribution 
fails to reproduce the data. 

Before addressing a risk analysis in the next Section, it is worth mentioning, for 
completeness, that other approaches to model with accuracy the tail exponent 
of the returns cdfs are discussed in the literature. They are based on Extreme 
Value Theory [18] and Hill's estimator [19,20]. However, since they mainly 
focus on the tails, they require very long time series to accumulate sufficient 
statistics and are not considered in the present study. 
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Fig. 2. From top left clockwise: Autostrade SpA, Telecom Italia (from May 15 th 
2001 to May 5 th 2005), Mibtel and Mib30 (from March 27 th 2002 to March 13 th 
2005) P> of negative daily returns. Points represent historical complementary cdf, 
while dashed and solid lines correspond to Gaussian and Student fits, respectively. 
The parameters values of the fitted curves are detailed in Table 2. 



a, 






-R -R 

Table 2 

Mean m, volatility a, and tail exponent for Autostrade SpA, Telecom Italia, 
Mibtel and Mib30 time series, m and a are estimated from empirical moments, 
while v is obtained through a negative log- likelihood minimization as in Eq. (10). 



m 



a 



Autostrade 0.12% 1.38% 2.91 

Telecom -0.02% 2.23% 3.14 

Mibtel 0.02% 1.03% 3.35 

Mib30 0.02% 1.16% 3.22 



4 Risk analysis 



In this Section we present a comparison of the results obtained estimating the 
market risk throught VaR and ES according to different methodologies. The 
standard approach is based on the normality assumption for the distribution 
of the returns. For this case we are provided of closed-form solutions, Eqs. (3) 
and (4), that depend on the two parameters m and a. For the time series 
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under consideration, the effect of the mean, as shown before, is negligible, and 
the surviving parameter is the volatility a. Several techniques are discussed 
in the literature to model and forecast volatility, based on stochastic volatilty 
approaches [21], GARCH-like [22] and multifractal models [23]. They usually 
require very long time series (tipically 300 high frequency returns per day over 
~ 5 - 10 years) and are quite demanding from a computational point of view. 
As discussed in Section 3, we limit our analysis to 1000 daily data and we 
estimate the volatility using the empirical second moment. In order to avoid 
the problem of a uniform weight for the returns, RiskMetrics introduces the 
use of an exponential weighted moving average of squared returns according 
to the formula [12] 

1 - A N ~ x 

where A G (0, 1] is a decay factor. The choice of A depends on the time hori- 
zon and, for At = 1 day, A = 0.94 is the usually adopted value [12]. a t +i\t 
represents volatility estimate at time t conditional on the realized R. If one 
considers Eq. (11) as the defining equation for an autoregressive process fol- 
lowed by <J t+1 \ t (coupled with R t = a t e t with e t ~ i.i.d. (0, 1 )) , Refs. [24,25] 
provide reasons for the claimed good success of the RiskMetrics methodology. 
In order to relax standard assumption about the return pdf without loosing 
the advantages coming from a closed-form expression, we presented in Sec- 
tion 2 generalized formulae for VaR and ES based on a Student-t modeling 
of price returns. In this framework, the tail index v emerges as a third rel- 
evant parameter, which is possible to constrain using a maximum likelihood 
technique, as previously described. As a benchmark of all our results, we also 
quote VaR and ES estimates following a historical approach, which is a pro- 
cedure widely used in the practice. According to this approach, after ordering 
the N data in increasing order, we consider the [NV*] th return R^nv*]) as 
an estimate for VaR and the empirical mean over first [NT*] returns as an 
estimate for ES 1 . 



At a variance with respect to previous investigations [11,24], we also provide 
68% confidence level (CL) intervals associated to the parameters. In this way 
we can estimate VaR and ES dispersion. To this extent, we implement a boot- 
strap technique [26]. Given the N measured returns, we generate M = 1000 
sinthetic copies of R, {R^}, with j = 1, . . . , M, by random sampling with 
replacement according to the probability p = (1/iV, . . . , 1/N). For each R* we 
estimate the quantities of interest and we obtain bootstrap central values and 



1 The symbol [ ] stands for integer part, while Rrj\ is the standard notation for the 
j th term of the order statistic of R. Since iV > 1 we neglect the fact that the p th 
entry is a biased estimator of the p/TV-quantile, i.e. E[i?( p )] = p/(N + 1). 
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Table 3 

Parameters values and bootstrap estimates for the 68% CL intervals for the time 
series as in Table 2, 





m 


a 


°"t+i|t 




#(10) 


Autostrade 
Telecom 
Mibtel 
Mib30 


o-i2±8:8S% 

-0.02l°;$% 
0-02±g;g% 


i 00+O.O80/ 

2.23±g;ii% 

i-03±8:S2% 

1.161S;«% 


1.831°,;!*% 

o.72i8:i% 


o 14+O.2I 

o qrr+0.18 
°-°°-0.19 

o 99+0.15 
°- zz -0.16 


-3.5iij} : g% 

-6.14l?;l% 
-2.961°;^% 
-3.3318;™% 



confidence levels. For example, we use for the mean the relations 



i M I AT-1 

™Z = ^I>; with m^-g^ (12) 

and we define the 1 — 2a CL interval as [to* , TO*_ a ], with to* such that P(to* < 
to*) = a and a = a, 1 — a. For 68% CL, a = 16%. In Fig. 3 and Tables 3, 4, 5 
we quote results according torn- (to| — m*) + (to*_ q — to*,). In this way, we 
use the bootstrap approach in order to estimate the dispersion of the mean 
around the measured value, m. Analogously for all the other parameters. 

Table 3 shows central values and estimated 68% CL intervals for the daily 
returns series under study. These numerical results come from a straightfor- 
ward application of the resampling technique. It is worth mentioning that it is 
possible, and sometimes necessary, to use improved versions of the bootstrap. 
As a rule of thumb, we consider the boostrap approach accurate when, given a 
generic parameter, the difference between its empirical value and the boostrap 
central value estimate is close to zero and 68% CL interval is symmetric to a 
good approximation. In our numerical simulation, we measured a systematic 
non zero bias for a t+1 \ t and from Table 3 it is quite evident the asymmetry of 
R([nv*]) intervals for both Autostrade and Telecom data. We can, therefore, 
consider the corresponding CL intervals as a first approximation of the right 
ones, since bias and skewness corrections would require sofisticated and ad-hoc 
techniques [26], which are beyond the scope of the present work. 

In Fig. 3 we show VaR and ES central values and 68% CL bars for Autostrade 
SpA and Mib30, corresponding to 1% and 5% significance level and according 
to the four methodologies previously described. In Tables 4 and 5 we de- 
tail all the numerical results, including also Telecom Italia and Mibtel data. 
As already noted in Ref. [24], at 5% significance level Student-t and Normal 
approaches are substantially equivalent, but here such a statement sounds 
more statistically robust, thanks to the bootstrap 68% confidence levels and 
to the comparison with the historical simulation. At this significance level, 
we register for VaR a different behaviour between single assets and indexes. 
While assets show the best agreement between the Student-t and historical 
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Fig. 3. VaR A* (upper panel) and ES E* (lower panel) central values with 68% 
CL intervals for Autostrade SpA (left) and for Mib30 (right), according to the four 
different methodologies discussed in the text. The significance level V* is fixed to 
1% (circles, solid lines) and 5% (triangles, dashed lines). 
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Autostrade Mib30 
approaches (see also Table 4), for Mib30 and Mibtel data we observe the best 

agreement between the Normal and historical methodology. In order to en- 
force this empirical evidence, it would be necessary to analyze additional time 
series to see to what extent this difference between assets and indexes holds. 
From Fig. 3, Table 4 and Table 5 it can also be seen that A* and E* central 
values calculated according to RiskMetrics methodology are quite fluctuating 
and characterized by the largest CL bars. The decreasing of V* traduces in 
a major differentiation of the different approaches. In general, we obtain the 
best agreement between the Student-t approach and the historical simulation, 
both for A* and E*, whereas, as before, the RiskMetrics methodology overesti- 
mates or underestimates the results of the historical evaluation and is affected 
by rather large uncertainties. 



To conclude, we would like to note that we expect, from the results shown in 
Fig. 1 and Table 1, that, for a fixed significance level, there exists a crossover 
value, Across, below which the generalized Student-t VaR and ES formulae un- 
derestimate the Gaussian predictions. This effect was already mentioned in 
Ref. [11], but the analytical formulae here derived allow us to better charac- 
terize it. Under the hypothesis of a Student-t distribution, the crossover value 
does not depend on the first and second moments and, therefore, the knowl- 
edge, for a given time series, of the tail exponent only is sufficient to conclude, 
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Table 4 

Estimated VaR values (mean and 68% CL interval) for 1% and 5% significance 
levels from Autostrade SpA, Telecom Italia, Mib30 and Mibtel. For each time se- 
ries, the results of Student-* and Normal fit, historical simulation and RiskMetrics 
methodology are shown. 





Student-t 


Normal 


Historical 


RiskMetrics 


Autostrade VaR 1% 
VaR 5% 


3.47210:1! 

1 71 7+0.071 
i -' i '-0.071 


o n qi+°- 197 
o.uyi_ 204 

9 1 sn+0.139 
Z.lOU_ .145 


3 516 +a149 
O.dlU_ 306 

i-8io±8:iH 


4.13818;^ 

2.890l8;il8 


Telecom VaR 1% 
VaR 5% 


5.9001°;™ 

o 191 +0.135 

°- izi -0.137 


5. 2001^277 

q fifi9 +0.214 
O.OOZ_ 202 


fi ,07+1.348 
u - 10 ' -0.866 

o oqo+0.110 
o.o»o_q 127 


q CQC+0.990 
°- oyo -1.085 

2.548l8; 6 7 ? 4 7 


Mib30 VaR 1% 
VaR 5% 


n 47 +0.106 
J - U4 ' -0.105 

1 «1 9+0-066 


2 fi7^+ - 097 
1 S8^+ - 073 


q qqi +0.255 

°- ool -0.304 

2.01018;™? 


1 662+ - 524 
r-ooz_o 516 

1 lfiO+0.375 
l.lDy_03 58 


Mibtel VaR 1% 
VaR 5% 


9 71 s+0.097 
z. i lo_ 092 

i 4 c 4 +0.062 

i.+o+_ 062 


q 7Q +0.088 
z "J'°-0.084 

1 fi74 +0.065 


2.96718;!° 
1.81ll8;f 7 ° 


1 +0-453 
l.O51_ .449 

1 iio +0 - 324 
i.iiu_ 316 



Table 5 

Estimated ES values (mean and 68% CL interval) for 1% and 5% significance levels. 
Time series and methodologies as in Table 4. 





Student-t 


Normal 


Historical 


RiskMetrics 


Autostrade ES 1% 
ES 5% 


5.503l8;4 2 } 

2.94618;11 


q cc-Q+0.229 
O.Ody_ .231 

9 797+O.I75 
Z -' Z( -0.182 


5.07618;^ 
3.006l8;ll 


4-75918J3 7 

q fic;c:+0.653 
O.OOO_ .677 


Telecom ES 1% 
ES 5% 


8.912l8:I 8 l 

K noc+0.242 
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9.68511; 4 ! 

r q 9n +0.478 
d.ozu_o.466 


4 1 1 fi++ 133 
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3 iqn+ - 879 

O.iyu_ .969 


Mib30 ES 1% 
ES 5% 


a c-79+0.199 
+- ' z -0.191 

9 5OR+0.093 
z.o»u_ og! 
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a priori, whether the fat-tailed results for VaR and ES will underestimate or 
not the corresponding Gaussian estimates. 



5 Conclusions 

In this paper we have presented a careful analysis of financial market risk 
measures in terms of a non-Gaussian (Student) model for price fluctuations. 
We have derived closed-form parametric formulae for Value at Risk and Ex- 
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pected Shortfall that generalize standard expressions known in the literature 
under the normality assumption and can be used to obtain reliable estimates 
for the risk associated to a single asset or a portfolio of assets. The obtained 
non-Gaussian parametric formulae have been shown to be able to capture ac- 
curately the fat-tailed nature of financial data and, when specified in terms of 
the model parameters optimized by means of an empirical analysis of real daily 
returns series, have been found to be in good agreement with a full historical 
evaluation. Moreover, the risk measures obtained through our model show non 
negligible differences with respect to the widely used Normal and RiskMetrics 
methodologies, indicating that the approach may have helpful implications for 
practical applications in the field of financial risk management. We also pro- 
posed a bootstrap-based technique to estimate the size of the errors affecting 
the risk measures derived through the different procedures, in order to give a 
sound statistical meaning to our comparative analysis. 

As far as possible perspectives are concerned, it would be interesting to in- 
vestigate to what extent our conclusions, drawn from an analysis of a sample 
of Italian financial data, apply also to other financial markets. In particular, 
one could check whether, at a given significance level, statistically relevant 
differences are present between the results valid for a single asset and those 
relative to a portfolio of assets, as our analysis seems to indicate, at least for 
5% VaR. Another interesting development concerns the comparison between 
the predicitions for VaR and ES of our model with the corresponding ones 
derived by means of other statistical procedures to measure tail exponents 
known in the literature [18,19,20], as well as with the results from simulations 
of advanced models of the financial market dynamics, such as GARCH-like 
and multifractal models [22,23]. 
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